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ABSTRACT 


Mathematical relationships between subjective comfort and 
environmental variables in an air transportation system are 
investigated. As a first step in model building, only the motion 
variables are incorporated and sensitivities are obtained using 
stepwise multiple regression analysis. The data for these models 
have been collected from commercial passenger flights. 

Two models are considered. In the first, subjective comfort 
is assumed to depend on rms values of the s ix-degrees -of-f reedom 
accelerations. Variations in the model for different subjects 
are also explored. 

The second assumes a Rus tenburg-type human response function 
in obtaining frequency-weighted rms accelerations, which are used 
in a linear model. The form of the human response function is 
examined and the results yield a human response weighting function 
for different degrees of freedom. In addition, an improved 
subjective comfort model is obtained. 



Section 1 


INTRODUCTION 

With the advent of short-haul aircraft over the last few years, 
increasing attention has been directed at aircraft ride quality inasmuch 
as many of these vehicles fly at low altitudes and encounter turbulence 
over a significant portion of their flights. The objective of this 
study is to develop a subjective comfort model which will predict human 
reaction to the s ix-degrees-of — freedom of aircraft motion, viz., three 
linear accelerations (vertical, transverse, and longitudinal) and three 
angular accelerations (pitch, roll, and yaw). The data utilized in 
obtaining the models contained herein are from the first phase of a 
continuing study. As future data is obtained refinements will be 
reported. 

1 . 1 Background to the Problem 

Host of the work in ride quality has been done since the 1930's — 
primarily in the field of automobile and railway ride comfort— and has 
been operator-performance oriented, rather than passenger-comfort oriented. 
The work that has been done in developing comfort criteria is summarized 
by Jacobson (1). Some of the vertical acceleration criteria are shown 
in Figure 1.1 and similar results for transverse acceleration are shown 
i n Figure 1.2. 

The major portion of the experiments done to date have used ground- 
based simulators although increasing attention has been directed at actual 
field testing (e.g., Jacobson (3) or Koo (4) ) in the last decade. Simulators 
have certain advantages; one of the most important is the ability to collect 
considerable data using a number of test subjects with minimal cost. Further, 
there is the advantage of wider control over the acceleration input (all the 
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FIGURE 1.1 VERTICAL VIBRATION CRITERIA 
(Taken from Jacobson (1) and Draft International Standard (2)) 
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FIGURE 1.2 TRANSVERSE VIBRATION CRITERIA 
(Taken from Jacobson (1) and Draft International Standard (2)) 
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s j x _degrees-of-f reedom) which is essential to obtain a good mathematical 
model. On the other hand, ground-based simulator tests lack the ability 
to simulate psychological factors (e.g. anxiety and motivation) among 
test subjects. This is especially undesirable for aircraft comfort 
tests, since the anxiety factor of being airborne cannot be simulated on 
the ground. 

Further, most previous work has involved human response to vertical 
acceleration only. Some recent research has been directed at other axes 
of acceleration such as transverse and longitudinal (e.g. ISO (2) and 
Koo (4)), but surprisingly little has been done on angular accelerations 
(e.g. Collins (5) and Clark and Stewart (6)). 

In addition, most of the simulator tests have been conducted with 
sinusoidaloscil lat ions at different discrete frequencies. Test subjects 
are generally exposed to increasing acceleration at each frequency and 
they rate the acceleration based on a comfort scale. The exposure limits 
given by the criteria are probably applicable only for pure sinusoidal 
acceleration in the corresponding degree of freedom. The results may not 
be applicable, if there is more than one frequency of acceleration occur- 
ring simultaneously (some efforts have been made in this direction-- 
Brumaghim (7) and O'Massey, et al. (8)) t or if there are multiple degrees 
of freedom present (an area where hardly any work has been done at all). 

Though most of the experiments were conducted under different physical 
environments (e.g. sitting, standing, etc.), the results indicate that 
humans are most sensitive to the frequency range of 4 - 8 Hz in the vertical 
mode, probably due to internal organ resonance. Except on this point, 
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there is wide disagreement on the magnitudes of acceleration for a given 
semantic level (e.g. perceptible, uncomfortable, etc.). The intolerable 
level of one investigator is not even mildly annoying for another (l). 
Such discrepancies are found widely, and experimental and semantic 
differences do not always explain the spread in results. 

Table- 1.1 indicates the sources of vibration with dominant 
frequencies, flight phases, and durations for typical aircraft. Of the 
various flight phases, cruise is the longest and the main sources of 
vibration are the engine (propeller or jet), turbulence, and buffeting. 

1 . 2 Statement of the Problem 

Vibration has been found to be an important variable affecting 
passenger comfort (3). Figure 1.3 shows the frequency ranges of 
mechanical vibration effects on man. It is seen that most of the effect 
on humans is caused by infrasonic frequency vibrations (0 to 25 Hz), 
hence this study is restricted to this range, 

Data was collected on regular commercial flights, with test subjects 
riding as passengers, creating both the proper physical and psychological 
environment as well as providing representative accelerations encountered 
i n regular f 1 ights . 

Several models are evaluated: 

a. A comfort model depending on the root mean square (rms) values 
of the six measured accelerations is developed. Inter-aircraft and 
inter-subject differences are explored. 

b. A comfort model depending on f requency-weighted rms accelera- 
tions using a Rustenburg-type human response function is developed. 
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TABLE 1*1 


AIRCRAFT MISS' 

ION DESCRIPTION AND 

PROMINENT VIBRATION SOURCES 


(Taken from O'Massey (8)) 


Flight Phase 

Approximate 
Time Duration 

Prominent 
Vibration Sources 

Frequency 

Range 

Warm up 

1 to 15 min. 

Engine 

40 and above 

Taxi 

5 to TO min. 

Runway Roughness 
and Engine 

0.5 to 5 
40 and above 

Run up 

2 to 15 min. 

Engine 

40 and above 

Takeoff 

1 to 3 min. 

•>v 

Runway, Turbulence 

Engine 

Buffet 

0 to 10 

40 and above 

1 to 20 

Climb 

3 to 30 min. 

Engine 

Turbulence 1 

40 and above 
0 to 10 

Cruise 

15 min. to 8 hrs. 

Engine 

Turbulence" 

Buffet 

40 and above* 

0 to 10 

1 to 20 

Descent 

5 to 15 min. 

■,v 

Turbulence 

0 to ?0 

Landing gear 
and flaps down 

1 to 10 min. 

Flap buffet and 
Turbulence" 

0 to 10 

Landing 

5 sec. to 2 min. 

Landing impact 
Runway Roughness 

Pulse in nature 

0.5 to 5 


Here turbulence refers to both atmospheric and boundary layer turbulence. 


6 


0.1 

- J — 


FREQUENCY 

10 1 


10 * 


- Hz 
. J 


10 


10 ^ 


10 ! 


I nf ra son i c 
vi brat 1 oris 


HUMAN RESPONSE 


Range of vibratory sensation 
Disequilibrium (high intensities) 
Motion sickness 
Major body resonances 
Blurring of vi sion 
Disturbance of speech 
Interference with manual tasks 
Discomfor t/pal n (high levels) 
General s tress/"f a t I gue" effects • 
Damage risk to spine* 
Gastrointestinal disorders* 
Vibration syndrome (hands)* 


* Cumulative occupational hazard 


Audio-frequency 

range 


Ultra- 

sonics 


■♦-SCOPE Of 
ISO 

LIMITS 


uncertain ranges 


FIGURE 1.3 APPROXIMATE FREQUENCY RANGES OF PRINCIPAL 
EFFECTS OF MECHANICAL VIBRATION ON MAN 

(Taken from Broderson and Von Gierke (9)) 
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c. A comfort 'model .depending, on frequency-weighted rms accelerations 
is determined. This model yields human sensitivity as a function of 
frequency for each degree of. freedom. 
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Sect ion 11 

DATA ACQUISITION AND REDUCTION 

2 . 1 Introduction 

A theoretical approach to the prediction of subjective comfort 
level, given, s i x-degrees -of-f reedom of motion, has been ruled out 
because of the lack of understanding of the complex interaction of 
biodynamic, physiological, and psychological characteristics of humans. 
The only recourse is to employ empirical techniques which require a 
large amount of data collection. 

Past ride quality studies have used a variety of comfort scales 
(1), most having 3 to 7 levels of comfort. Since people prefer not 
to use end points , the 3-level scale is inadequate. For the present 
work, a 5-level comfort scale has been adopted : 

1. Very comfortable 

2. Comfortable 

3. Neutral 

A. Uncomfortable 
5. Very uncomfortable 

No performance measures have been associated with these levels since 
test subjects might be misled in -evaluating performance measures rather 
than comfort. 

2.2 On-Board Data Recording 

The measuring and recording system, designed and fabricated by 
NASA Langley Research Center (shown in Figure 2.1), consists of battery 
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FIGURE 2.1 PORTABLE INSTRUMENT PACKAGE 
AND RECORDING EQUIPMENT 
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operated sensors and recording equipment weighs a total of 30 pounds. 
Linear accelerations to + 1/2 g in each of three mutual ly-perpendicular 
directions (vertical, longitudinal, and transverse), angular accelerations 
to + 2 rad/sec about three mutual ly-perpendicular axes (pitch, roll, 
and yaw), subjective response, and reference frequency signal are all 
recorded on a two-channel 1/4" magnetic tape recorder (AH) --one channel 
of FH multiplexed data, the other voice. A functional block diagram 
of the measuring and recording system is shown in Figure 2.2. The 
voltage-controlled oscillator (VGO) for each of the recorded variables 
corresponds to standard inter-range instrumentation group (IRIG) 
frequencies (10), (11), and (12). The reference frequency signal is . 
used during data reduction to compensate for tape speed variations. 

The accelerometers are placed on the floor of the vehicle, at the 
foot of the subject's seat. The subject indicates his comfort level by 
depressing the appropriate comfort box button, corresponding to one of 
the semantic levels described above, when evaluating the ride over a 
period of approximately 30 seconds, in 2-minute intervals during flight. 

2.3 Data Recovery 

A functional block diagram of the data recovery system is shown 
in Figure 2.3. The multiplexed data processed through discriminators 
has an output which is directly proportional to the input signal 
frequency and hence to the amp! i tude of acceleration. 

During playback, the reference frequency is used by a compensation 
module which detects deviations due to tape speed changes and applies an 
inverse percentage compensation, considerably reducing errors. Each. 
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FIGURE 2.2 FUNCTIONAL BLOCK DIAGRAM OF MEASURING 
AND RECORDING SYSTEM 
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FIGURE 2.3 FUNCTIONAL BLOCK DIAGRAM OF DATA 
REDUCTION SYSTEM 
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channel of data is recorded on an eight-channel wide-band FM tape with 
carrier frequency (6.75 KHz) along with a continuous time code. This 
time code becomes the primary reference in the data analysis. An 
oscillograph record of all measurements complete with the time code, 
subjective responses, and acceleration signals is obtained simultaneously. 
A typical output, less the time code, is shown in Figure 2.4. The 
recovered data is processed through a digitizer, sampled every 0.02 
second (each channel), and written on a standard digital tape. 

2.4 Data Reduct ion 

It is nearly impossible to use acceleration time histories alone 
to arrive at a subjective response model. The data must be reduced to 
meaningful variables which are amenable to modelling. A Time Series 
Analysis program (TSA) is utilized at the LRC computer facility for 
providing various statistical quantities (13). These include the 
means, standard deviations, autocovar iances , crosscovar lances , and 
power spectral density functions for each of the s Ix-degrees-of - 
freedom. These quantities, along with a ten-digit code consisting of 
the date of flight, type of aircraft, flight number, and a sequence 
number uniquely identifying a given flight segment, are stored on 
magnetic tape to be used in the data analysis. 
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Section III 


ANALYSIS OF DATA 

3 . 1 Assumptions In Developing the Model 

Most of the vibration experienced by the test subject is through 
the seat, thus there is a discrepancy with the vibration measured by 
the equipment at the foot of the test subject. For the present study, 
it is assumed that the seat transfer function is the same for ail the 
aircraft under cons iderat ion in the frequency range of interest. 

3.2 Data Col lection 

All the data used in this study were taken on commercial airline 
flights between January and April of 1972. This involved approximately 
100 flight segments flown on-board three different aircraft — F-227, 

YS -11, B~737- ta under a variety of turbulence conditions and over 
dissimilar terrain. (For aircraft details see Appendix A.) The 
number of flight segments per day varied from 2 to 8, from approximately 
60 miles to 250 miles in length, and from approximately 15 minutes to an 
hour in duration. One or two test subjects flew on-board each flight 
segment and a total of nine test subjects participated in the program. 
(See Appendix A for details.) 

3.3 Mode 1 1 i ng 

Regression analysis is used to develop subjective comfort models 
(a detailed description of the statistical quantities used, as well as 
the regression analysis can be found in reference 14) where subjective 
comfort is the dependent variable and accelerations in the six-degrees- 
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of -freedom are the Independent variables. The regression program 
used is part of the UVA system library (SPSS, Statistical Package 
for the Social Sciences) (16). 

The first model developed uses autocovar I ance and crosscovariance 
and has the form: 


C = b. 


6 

+ I b 

J=l 


J (S j 5 


6 6 

+ I I 

j=l k-j+l 


b jk {i j k ) 


(3.0 


where 


C - predicted value of the subjective comfort 


b-,b., - constants predicted by the regression analysis 

J j k 

a. = /r v (0) - rms value of the acceleration 


J 


(3.2) 


J 


a., - j/r y y (0 ) - zero lag crosscovariance of variables X. and X, • (3.3) 

J k a A Ik 


J k . 


Hence, from Equations (2.5) and (2.7): 


a/ = j / X 2 ( + ) dt, 
0 J 

and 


(3.4) 


■ T 

^j k = 7 / Xj(t) X k (t) dt. (3.5) 

With appropriate modification nonlinear terms such as square roots can 
also be incorporated in the above model. 

The second type of model uses a Rustenburg-type (17) human response 
function to compute frequency-weighted rms accelerations. It is necessary 
to obtain a human equal sensation curve; this is not easy, however,., 
because existing data have large scatter as is seen in Figures 1.1 and 
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and 1.2. Since most of the available data is for subjects exposed to 
vertical vibration, a human equal sensation curve has been adopted based 
on that data and is shown in Figure 3.1. (Taken from Jacobson (1).) 

The human frequency response function, W(f ) , is obtained by inverting 
the equal sensation curve as is shown in Figure 3.2. W(f) indicates the 
relative weight at frequency f that subjects use in evaluating comfort, 
and is normalized over the entire range of frequencies of interest. 

An assumption made in this model is that the response function 
W(f) is valid in all s ix-degrees-of -f reedom, This assumption is 
inherently weak as can be seen in Figures 1.1 and 1.2, but is used in 
the hope that it will indicate trends. 

The frequency-weighted rms acceleration (of X.) is given by: 


where 


= / W 2 (f) ^ . (f) df 

j o j 


F - the upper limit of frequency. 


( 3 . 6 ) 


To normalize W(f), ^ (f) is set equal to 1 in Equation (3.6) and 

j 

the weighted and unweighted (W(f) - 1.) cases are compared, yielding: 


'/ W 2 (f) df * F. 


( 3 . 7 ) 


W (f ) is gotten by interpolation as: 

wo = w (y (^y' k 


( 3 . 8 ) 


for f k 1 f 1 f k+ l 


where ^ is the slope between the break points of the W.(f) curve. 
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Table 3.1 tabulates the values of ^ used for the model ; 


C = b + T b. (a. ) 

° j=i J w j 





TABLE 3.1 

SLOPES 

OF THE 

HUMAN RESPONSE FUNCTION, W (f) 

Ranqe of 

Frequencies (Hz) Slope ^ 

0 

to 

1.5 

0 

1.5 

to 

4 

1 

4 

to 

8 

0 

8 

to 30 

-1 


(3.9) 


The final model considered for analysis is frequency dependent, 
and has the form: 


6 K 


C = b 


o' * X I b 'ik <5 i V 
j=l k= I Jk J ’ k 


( 3 . 10 ) 


where 

C - predicted value of comfort 
b’, b’ ( - coefficients predicted by the regression analysis, 

O J K 

and where the rms value is given by (15): 

f, 


k+l 

5 j,k = { V <f) df ' 

f k J 


(3.11) 


Here 


4>v (f) - p.s.d.f. of the variable X. 


f - frequency 
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f^ + | - lower and upper limit of the frequency band, 

over which the rms a. . is evaluated 

J» k 

K - number of frequency bands used. 

Suitable frequency bands are chosen in order to develop this model, 
yielding a frequency-dependent model and a human response function for 
each degree of freedom. 

These three models serve as the preliminary criteria of ride 
qua 1 i ty. 
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Section IV 

RESULTS AND DISCUSSION 


4. I The linear Model 

Models are restricted to a few coeff ici ent$--about 100 data segments 
being needed for each. By this criterion, the largest model can have 14 
coefficients at most (approximately 1400 data segments are available). 

As a first step, only the rms values derived from the autocovariances 
are used in a linear model, resulting in a comfort response, C, which 
is related to the six accelerations by:* 

C = 1.82 + I4..3a v + 0.l63a R + t.0a L + 2.3a T + 0.2a y + 0.l4a p . (4.1) 


Vertical is found to be the dominant and best defined variable 
for the following reasons: 

a. Vertical acceleration accounts for the largest average contri- 
bution to the comfort value compared to the other degrees of freedom. 
The average contribution of each acceleration component, Xj , is 
computed as follows: 


where 




(4.2) 


C- is the mean contribution to C by a. 

a j J 

y- Is the mean value of a., 
a . 1 


The average values of C- are: 

j 


’•In this and ensuing models, all linear accelerations have units of rms 
g's and all angular accel erat ions rms rad/sec 2 . 
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Vertical 0.92 

Longitudinal 0.073 

Transverse 0.048 

Roll 0.051 

Pitch 0.022 

Yaw 0.011 


The comfort scale is divided into five step function values (i.e., 
integer values of 1 to 5 ) implying that any degree of freedom which 
on an average contributes less than 0.5 stands less than 50 % chance 
of altering the step response "C" by itself.* Hence all such variables 
are relatively insignificant individually (in comparison with those 
that contribute more than 0.5 on an average). In this sense only 
vertical seems to be significant. 

b. As can be seen from Table 4.1, the correlation of vertical 
acceleration with the comfort rating (0.723) is the largest. 

c. Further it is seen that the coefficient of variability (which 
is defined as the ratio of the standard error to that of the expected 
value or mean) of the predicted coefficient of vertical is low. The 
variability coefficients are listed below: 


Vertical 

0.04 

Roll 

0.51 

Long i tud i nal 

1.72 

Transverse 

0.588 


-Here the assumption Is made that human response varies continuously 
between discrete levels of C and that the crossover occurs at the 
midpoint. 
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TABLE A. 1 


Comfort 

P i tch 

Rol 1 

CORRELATION MATRIX (SYMMETRIC) 

Yaw Trans. . Long. Vert. 

v Vert . 

T rans .- Vert . 

Comfort 1.0 

0.107 

0.209 

0.273 

0.53 

0.271 

0,723 

0.724 

0.544 

P i tch . 

1.0 

0.836 

0.366 

0.2 

0. 123 

0.031 

0.018 

0,12 

Roll 


1.0 

0.442 

0,351 

0 . 18 1 

0. 164 

0.15 

0.-207 

Yaw 



1.0 

0.472 

0.197 

0.342 

0.315 

0.354 

Transverse 




1.0 

0.351 

0.751 

0.726 

0.699 

Long 1 tudinal 





1.0 

0.437 

0.454 

0.369 

Vertical 






1.0 

0.982 

0.763 

Alert ical 







1.0 

0.729 


Trans.- Vert. 


1*0 



Yaw 


9.59 


Pi tch O .96 

The lower the coefficient of var tab i 1 1 ty, the greater the confidence 
in the coefficient. Hence, the greatest confidence can be placed on 
the coefficient of vertical. 

d. The model, Equation (4.1), can be expressed in a normalized 
form as fol lows : 

C 1 = 0.71a v ' + 0 . 07 a R ' + 0 . 02 ^' + 0.05^' + 0 . 02 ^’ + 0.035a p ' 

(4.3) 

where 

C = (c - c)/a c 

C - mean of C 

a' - rms values based on autocorrelat ion with 

zero mean and unity standard deviation over 
the set of n data cases 

T - standard deviation of C. 

W ‘i 

Here each of the independent variables have zero mean and unity 
standard deviation. If It is assumed that all of the normalized 
variables have a similar probability density distribution, then the 
probabilities of each variable taking a given value are almost equal. 

Since the coefficient for a ^ 1 is the largest, the average contribution 
from a^’ to C', and in turn to C, is the largest, again Indicating that 
vertical is the most dominant. 

e. Finally, the r value (see reference 15 for definition) indicates 

that most of the improvement was achieved with vertical alone in the model 

7 2 

(r = 0.52) versus addition of all other terms (r = 0.54). 
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4.2 The Nonlinear Model 

Since vertical is the dominant variable in determining subjective 
comfort, it is given an additional degree of freedom. C - 5 represents 
the highest comfort level independent of the acceleration. Thus a model 
which linearly increases for small accelerations and flattens out at 
high accelerations would be preferable. Incorporation of the additional 
term v/3^", modifies the linear model to behave in the desired fashion, 
yielding! 


C = 1.33 + 5.29 ^ + 0.158a R + 1,0a + 4,0i v + 2.69a T + 0.l6a 

+ 0. 09a Y • 


(4.4) 


The mean contributions to C in this case are: 


/Vert ical 

1.27 

Vertical 

0.26 

Transverse 

0.056 

Longitudinal 

0.073 

Rol 1 

0.05 

P i tch 

0.025 

Yaw 

0.005 


Once again vertical is the dominant variable. Comparing the 
above with the mean contributions of each variable in Equation (4.1), 
it can be observed that /I^* has replaced as the term contributing 
most. This is encouraging, since rises slower than and hence 
the model behaves better, for large values of vertical acceleration. 
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4.3 Model Using Vertical and Transverse Dependent Variables 

From Table 4.1, it Is seen that transverse has the second highest 
correlation with the comfort index. However, it is not the second 
largest contributor to C. This is due to the strong correlation between 
vertical and transverse (0.75), which is also confirmed by Figure 4.1 
showing the interdependence of vertical and transverse. The contribution 
of transverse to C occurs in two ways: partially through the transverse 

term which is independent of vertical and partially by vertical to which 
transverse is highly correlated. 

From Table 4.1 it is seen that only transverse, vertical, and rms 
crosscovariance between transverse and vertical (written for simplicity 
as trans-vert leal ) have a correlation wi th subject ive comfort of more 
than 0.5. Using these three variables, the model; 

C = 2.0 + 13. 8a y - 3.25a VT + 4.5a T (4,5) 

is obtained. Here note that the coefficient for is negative. This 
implies that the presence of each degree-of rfreedom masks the other to 
some extent. A similar effect is seen by Brumaghlm (7), in his study 
of dual frequencies. 

Depending on the application, one of Equations (4.1), (4.4), or 
(4.5) can be used for a criterion. Equation (4.5) is the simplest, 
whereas Equation (4.4) is the most accurate and complex model of the 
three. In most cases. Equation (4.5) should be adequate. 

4.4 Rustenburg-Type Model 

The second model is the Rustenburg-type (17) using frequency- 
weighted rms acceleration as shown in Equation (3.9) that yields: 
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TRANSVERSE 



LEGEND 


CCMf-l 

C0MF=2 




CC^F=3 
CGMF =m 



VERTICAL 


FIGURE 4.1 PLOT OF OBSERVATIONS ON VERT I CAL -TRANSVERSE PLANE 




(4.6) 


C = 1.31+ 6.54 vC7 + 0.46i + 0.5i w + 5.63a + 0.2 i 

W V R W L V Y 


+ 0,71a + 0.48a. . . 

W p w T 


The mean contributions to C of each variable is: 


/Vert tea i 

1.26 

Vert leal 

0.23 

Transverse 

0.0065 

Long i tudi nal 

0.005 

Pitch 

0.036 

Roll 

0.072 

Yaw 

0.008 


Equations (4.4) and (4.6) compare favorably except for the 
coefficient of transverse acceleration, which has decreased considerably. 
This is also confirmed by comparing the mean contributions of the two 
models. One of the reasons for the discrepancy between the two models 
is that the response function (or the weighting function), W(f) , is 
not a good description for transverse (this will be confirmed below). 

Table 4.2 shows the difference in correlation coefficients, p, 
between subjective comfort and the frequency unweighted and weighted 

A 

accelerations. As can be seen, the correlation between C and frequency- 
weighted rms accelerations is greater than- the correlation with the 
unweighted rms accelerations in all the degrees of freedom, indicating 
a definite frequency dependence of the comfort model. The same conclusion 
is reached by comparing the standard errors (0.568 for weighted and 0.573 
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TABLE 4.2 


COMPARISON OF CORRELATIONS, p, BETWEEN 

THE FREQUENCY 

UNWEIGHTED AND 

WEIGHTED VARIABLES AND 

COMFORT INDEX 


Independent Variable (X^) 

Unweighted 

p Ca. 

J 

Weighted 

p *w. 

J 

Pitch 

0. 107 

0.436 

Rol 1 

0.209 

0.333 

Yaw 

0.273 

0.355 

Transverse 

0.58 

0,593 

Longitudinal 

0.271 

0.356 

Vert ical 

0.723 

0.728 

/Vertical 

0.724 

0.728 
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unweighted) and the percentages of points with error greater than 0.5 
(39.6% for the weighted against 40.5% for the unweighted case). 

4. 5 Frequency-Dependent Model 

In order to obtain a better representation of human sensitivities 
than that given by Figure 3.2, a general frequency-dependent model 
(Equation (3.10)) is used. 

Only the p.s.d.f. in the range 0 to 8 Hz is used since in almost 
all instances, most of the energy is concentrated in this range. A 
three-band frequency split, 0 to 1.5 Hz, 1.5 to 4.0 Hz, and 4.0 to 8.0 
Hz, is chosen. 

If all degrees of freedom, each with three frequency splits, are 
incorporated simultaneously, the number of unknowns, 22, exceeds that 
which can be determined for the amount of existing data. Thus human 
sensitivities in each degree of freedom are considered individually. 

Vertical yields : 

C = 2.07 + 8.0a v , + 2 + 20; 3 ( 4 *7) 

where 

a. . - rms acceleration for 0 to 1.5 Hz 
” » ' 

a - rms acceleration for 1.5 to 4 Hz 

V , 4 

a.. - - rms acceleration for 4 to 8 Hz. 

The coefficients of ^ and ^ imply that subjects are most 
sensitive to frequencies in the 4 to 8 Hz range and least sensitive to 
frequencies in the 0 to 1.5 Hz range. This confirms the form of W(f) 
shown in Figure 3.2, and also agrees in relative magnitudes. 
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(4.8) 


For transverse: 

C = 2.4 + !5i Tj] + 47i T _ 2 + 4.7i Ti3 - 

As Is seen, a T (1.5 - 4 Hz) is the most heavily weighted whereas 
a T is the least weighted. This is in keeping with the data in 
Figure 1.2, and indicates that the response function, W(f ) , shown in 
Figure 3.2 does not hold for transverse acceleration. 

Other degrees of freedom are not evaluated due to lack of spread 
in the data, i.e. the mean acce Jerat ions in other degrees of freedom 
are nearly constant and hence their relationship with comfort cannot 
be established with any confidence, 

4.6 Inter-Aircraft and Inter-Subject Differences in the Node! 

The data was acquired on board three aircraft — F-227, YS-11, and 
B-737--using 9 test subjects. Here inter-aircraft and inter-subject 
differences in the model are evaluated. The simplest model, using 
rms values based on vertical, transverse and crosscovariance between 
vertical and transverse, is used for the analysis. 

Figure 4.2 shows inter-aircraft differences, where numbers in 
parentheses are the standard errors in the predicted coefficients. As 
is seen, the coefficients of vertical for all aircraft agree closely 
with the overall model. The coefficient of transverse for the overall 
model agrees well with that of the YS-11, reasonably with that of the 
B -73 7 » and poorly with that of the F-227. Note the fact that the 
standard error in prediction for the transverse coefficient in the 
case of the F-227 is high. The same behavior is observed for the 
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OVERALL: 


C = 2.0 + IJ.S**- 635 ^ - 3.25 (±K89) a VT + 4, 5 2^ +1 * 29) a T 


C * 2.95 


AIRCRAFT F-227: 


C = 2.395 + l4.l8 (±3Jl4) i v + 7.27 (ifi * l 4 ) a VT - 1! . 14 ( 47,J|6) 5 T 


c = 3.25 


AIRCRAFT YS-1 I: 


i ~~ 2.01 + ,3.83 (±0 ‘ 68) o v - 3.78 (i 2 -° 8 ) i VT + 4.4 5 ^ 3 ^ T 


C * 2.975 


AIRCRAFT B -73 7 : 


C = 1.9+ 14.29^'°^.- 4.04^ ±5 ’ 83 ^a VT + 10.15^*^a T 


C = 2.77 


Numbers in parenthesis indicate the standard error in predicted coefficient 


FIGURE 4.2 INTER-AIRCRAFT DIFFERENCES IN THE MODEL 
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coefficient of trans-vert ical . This discrepancy is probably due to the 
sparsity of data taken on-board the F-227. Further we see that for the 
flight test conducted the F-227 appears to be the most uncomfortable 
(C = 3.25) and the B-737 (C = 2.77) the most comfortable among the 
three aircraft. 

Similarly Figure 4.3 shows inter-subject differences for 3 subjects. 
The coefficients of vertical are similar with the exception of subject C. 
A large standard error resulted in the elimination of the coefficient 
of transverse for subject C. The overall model's coefficient of 
transverse compares favorably with that of subject B and reasonably 
with that of subject A. The coefficient of seems to be well behaved. 
The dispersion exhibited by the model for subject C is due to a sparsity 
of data for this subject. . 

Some deviations are observed for the coefficients of and a^, 
however they do not contribute significantly to the comfort model as 
confirmed in Figures 4.4 and 4.5. 

It can be concluded that inter-aircraft and inter-subject 
differences are small and hence the overall model can be used to 
represent all aircraft and subjects. 
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OVERALL: 


C = 2. 

C = 2. 

SUBJECT A : 

C — 1 . 

C » 2, 

SUBJECT B ; 

A 

C = 2. 
C = 2 

SUBJECT C ; 

C = 2 

C = 2 
"Numbers 

FIGURE 


0 + l3.8 (±0 - 635) i v - 3.25 (±K89) a VT + 4.52^ 


95 



.99 + I3.93 (±2 ^ 5) i v - 2.34 (±7J5) a VT 


!l.24 (±5 '* 57) i 


98 


09 + 12.5 ^ ±0 ' Sl 3 a y - 2.09 (± 2 * 3 ) a VT + 4.08^ ± 1 , 6 ) a T 


.97 


,02 +16.36^ l - 78 )i v .5.9, {±6 - 84) a VT 


■ 93 


in parenthesis indicate the standard error in predicted coefficient 


4.3. INTER-SUBJECT DIFFERENCES IN THE MODEL 
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LEGEND 
o OVERALL 
F AIRCRAFT F-227 
Y AIRCRAFT YS-il 
B AIRCRAFT B-737 


<F»B 


fb 

oY 


v'A-ay c vfAo vt c t'Aa| 


FIGURE 4.4 INTER-AIRCRAFT DIFFERENCES IN THE AVERAGE 
CONTRIBUTION TO COMFORT "C" 
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LEGEND 
© OVERALL 
A SUBJECT A 
B SUBJECT B 
C SUBJECT C 





Section V 


CONCLUSION 

5.1 Some Remarks 

Three types of models have been developed: models using rms 

accelerations; a Rustenburg-type frequency-weighted rms acceleration 
model; and a general frequency-dependent model. 

It was concluded that human sensitivities obtained for vertical 
acceleration, unlike transverse, follow the general form of the 
response function shown in Figure 3.2. Further, it was observed that 
the Rustenburg-type model was superior to the rms acceleration model; 
however the human response function Is not known for all directions, 
hence this type of model cannot be used to jts full potential. 

As a criterion, the simplest model developed, i.e., Equation (4.5), 
is suggested: 

C * 2.0 + 13.8i v - 3.25a VT + 4.52£ T (4.5) 

where inter-aircraft and inter-subject differences are minor. It Is 
cautioned that this is a prel iminary model which will be refined as 
future data becomes available. 

5.2 Appl icat ions 

Among the many applications for Equation (4.5) is the evaluation 
of existing aircraft ride quality under various flying configurations. 
Those found objectionable can be modified to conform to the passenger 
comfort criterion. The ride criterion can also be used in designing 
aircraft and flight control systems. 

Once a model has been selected, a suitable value of C can be chosen 
as the comfort boundary- -one between 3 and 4 would be a good choice (18).. 
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(5.1) 


Taking a value of 3.5, which yields 68% passenger satisfaction (18) 
Equation (4.5) yields: 

13.8a v - 3.25a VT + 4.52i T = 1.5. 


The crosscovariance is related to the rms accelerations in vertical 
and transverse by: 





(5.2) 


where P is the correlation coefficient for transverse with vertical. 
The value of P^j is obtained from Table 4.) ( p ^ = 0.751) yielding for 
Equation (5.1): 


1 3.8i v - 2.816 ^ ^ + 4.52i T = 1.5. (5.3) 

A. ** 

Equation (5. 3). along with isocomfort curves for C = 3 and C = 4 
are plotted in Figure 5.1. 

Figure 5.2 shows the superposition of the comfort criterion 
(C “ 3.5) on the actual data and indicates that more than 50% of the 
data lie within the acceptable range. It is also felt that the model 
may not be valid for large transverse accelerations simultaneous with 
small vertical accelerat ions , since negligible data was obtained in 
that range, thus the model represents an extrapolation in this region. 
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RMS TRANSVERSE ACCELERATION IN g's (aT) 




TRANSVERSE 



VERTICAL 


FIGURE 5.2 COMPARISON OF COMFORT CRITERIA WITH ACTUAL DATA 
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APPENDIX A 


SOME REMARKS ON THE AIRCRAFT AND SUBJECTS USED 
FOR DATA COLLECTION 

This appendix lists special characteristics and conditions under 
which the data were collected. 

Data were collected on three aircraft — FH-227B-200, YS-11A, and 
B - 73 7-201 (abbreviated in the text as F-227, YS — 11, and B — 73 7 » 

respectively). Table A.l indicates their characteristics and pertinent 

\ 

conditions under which the data were obtained. As can be seen, most of 
the data have been collected on board the YS - 1 1 . 

Table A. 2 shows the subject breakdown of the nine male test 
subjects used. Most of the data were collected using subjects A, B, 




TABLE 

A. 1 


differences 

'• 'i 

IN THE AIRCRAFT 

/ 

USED FOR DATA 

COLLECTION 

f; 

A! rcraf t 
Length (ft) 

FH-227B-200 

YS-11A 

B— 737-20 1 

84 

86 

100 

Wing span (ft) 

95 

105 

93 

Capacity 

(No. of passengers) 

44 

60 

90 

Max. wing loading 
(lbs/sq. ft.) 

60.3 

52.8 

115.2 

Eng i ne 

2 turbo-prop 

2 turbo-prop 

2 turbo-fan 

Normal terrain 
of f 1 ight 

Flat 

Mountainous 
and Flat 

Flat 

Average duration 
of f 1 ight (min.) 

17 - 55 

11 - 60 

20 - 52 

Altitude 
(1000 ft.) 

5 - M 

2 - 10 

6 - 26 

Stage length 
of flight (miles) 

72 - 240 

45 - 260 

85 - 380 

Cruise T.A.S. 
(M.P.H.) 

280 

280 

320 - 550 

Noise (db) 

. (Cp weighted) 

90 - 98 

80 - 100 

75 - 105 

Temperature 

(°F) 

70 - 82 

65-80 

65 - 78 

Seat comfort 

Least leg 
room 

Average leg 
room 

Most leg 
room 

No. of flight 
segments 

10 

7! 

14 



TABLE A. 2 

DIFFERENCES IN THE SUBJECTS USED: FOR DATA COLLECTION 


Sub i-ect 

Approximate 

Age 

(Years) 

Occupat i on 

Number of 
Flight 
Seqments 

A 

30 

Uni vers i ty 
teaching /eng i neer 

17 

B 

25 

Research Engineer 

82 

c 

25 

Research Engineer 

24 

D 

24 

Research Engineer 

13 

E 

40 

University 

teaching/engineer 

6 

F 

40 

University 
teachi ng/eng t neer 

4 

G 

■ 30 . 

University 

teachi ng/mat hemat ic i an 

2 

H 

52 

Uni vers ity 
teaching/eng i neer 

2 

1 

30 

Un ivers ity 
. teach mg/eng i neer 

1 
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